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ABSTRACT 


Theoretical  and  experimental  statistical  analysis  of  the  random  response  of  a 
continuous  bar  in  tension  is  presented*  Particular  attention  has  been  paid  to  the 
probability  distribution  of  the  strain  response  which,  for  a  linear  second-order 
system  under  gauss ian  excitation,  follows  a  Rayleigh  distribution.  However,  when 
the  excitation  level  of  the  clamped-clamped  continuous  bar  is  sufficiently  high  so 
that  the  tensile  strain  becomes  comparable  with  the  bending  strain,  then  the  strain 
crest  distribution  no  longer  follows  the  Rayleigh  prediction.  At  high  strain  levels 
the  distribution  of  positive  crests  as  well  as  maxima  is  greater  than  the  Rayleigh 
prediction  and  the  distribution  of  negative  crests  as  well  as  minima  is  less.  The 
distribution  of  positive  maxima  falls  below  the  positive  crest  distribution  as  the 
Q  of  the  system  decreases.  Similarly  the  distribution  of  negative  minima  falls 
below  the  negative  crest  distribution  as  the  Q  decreases. 
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Section  I •  Introduction 


In  many  investigations  of  response  of  mechanical  systems  to  random  excitation 
determination  of  mean  square  response  (displacement,  velocity  or  strain)  is  deemed 
sufficient  as  an  end  in  itself.  This  is  reasonable,  since  often  only  these  quanti¬ 
ties  are  measured  in  the  laboratory;  hence,  at  most,  we  need  the  probability  distri¬ 
bution  of  the  values  of  the  random  variable  at  a  given  time.  However,  when  the 
extremal  statistics  of  the  response  are  of  interest,  i.e.  ,  the  distribution  of  the 
maxima  or  minima  of  the  response,  more  statistical  information  is  needed,  preferably 
the  simultaneous  joint  probability  distribution  of  the  random  variable  and  its  first 
two  derivatives.  Such  considerations  arise  when  fatigue  theories  are  to  be  developed 
and/or  checked.  Fatigue  theories  will  not  be  discussed  further  here  but  they  will 
serve  as  the  motivation  for  the  following  theoretical  and  experimental  investigation. 
The  reader  is  referred  to  the  literature  for  such  discussions  (refs.  1,2,3,18,19). 

In  order  to  analyze  the  vibrations  of  an  actual  mechanical  system  it  usually 
is  approximated  by  an  ideal  system  such  as  a  plate  or  a  bar;  the  bar  often  is  deemed 
sufficient  for  illustrating  the  type  of  results  to  be  expected  for  more  complicated 
systems.  This  report  will  deal  with  the  vibration  of  a  bar  with  a  dynamical  non¬ 
linearity  in  the  form  of  a  tensile  stress;  for  the  sake  of  analysis  it  will  be  as¬ 
sumed  that  there  is  only  one  mode  of  the  bar  set  into  vibration.  This  procedure 
parallels  that  of  Lyon  (ref.  ^),  where  the  resulting  differential  equation  for  the 
displacement  of  the  fundamental  mode  is  reduced  to  a  hard-spring  oscillator  equation. 
For  experimental  and  theoretical  reasons  to  be  elaborated  on  later,  the  displacement 
response  will  be  equivalent-linearized  so  that  the  strain  extrema,  consisting  of 
bending  plus  a  tensile  component,  can  be  studied.  The  approximations  hold  so  long  as 
the  single  mode  model  is  valid;  it  must  be  remarked  that  at  higher  strain  levels 
(i.e.,  larger  displacement  amplitudes)  various  assumptions  break  down;  the  mode  shape 
changes  from  the  ^linear"  shape,  mode  coupling  cannot  be  neglected  and,  moreover, 
the  modal  damping  must  be  considered  nonlinear  even  if  linear  stress-strain  relations 
are  assumed  (see  for  example,  refs.  13,1^).. 

In  this  work,  the  so-called  crest  statistics  are  studied  in  detail  although 
theoretical  results  also  are  giv  :n  for  the  extremal  statistics.  In  addition,  ef¬ 
fects  of  the  strain  nonlinearity  on  the  average  number  of  zero  crossings  and  maxima 
per  second  are  considered  theoretically.  It  was  found  theoretically  and  experi¬ 
mentally  that  the  distribution  of  the  maxima  of  strain  as  well  as  the  positive 
strain  crests  lie  above  the  linear  (gaussian)  Rayleigh  values  for  the  higher  stress 
levels  studied.  On  the  other  hand,  the  negative  minima  and  negative  crests  fall 
below  the  Rayleigh  curve. 
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j  Section  II.  Statistical  Theory  of  Continuous  Bar 

I  A.  Displacement  of  a  bar  with  tension 

j 

I  Effects  of  a  small  stress  nonlinearity  on  the  extremal  statistics  of  a  simple 

J  supported  bar  forced  by  gauss ian  random  noise  have  been  studied  theoretically  by 

1  Lyon  (ref,  4),  A  parallel  analysis  for  a  bar  clamped  at  both  ends  but  extended  to 

i  somewhat  larger  strain  nonlinearities  is  here  reported. 


form; 


The  equation  of  motion  for  the  displacement  of  a  bar  with  tension  is  of  the 


y(x,t)  =  F(x,t) 


(ii-i) 


where  the  tensile  stress  T^  is; 


T,  = 


(11-2  ) 


f 


For  the  fundamental  mode  of  the  clamped  bar  with  tension  we  have,  approximately: 


y{x,t)=  y,(t)'l',(x) 


(II-3) 


with 


'|',(x)=  A  (coshX.x  -  cosXjx)  — 0.9825 (sinhXjX  —  sinX,x) 


] 


as  explained  in  section  III  of  this  report. 

Substitution  of  Eq,  (II-3)  into  Eq.  (II-l)  yields  the  following  equation  of 
motion  for  y^(t)  ; 


where 


( +  2a  +  ^yf )  y,  =  f  (t) 


dx)  =  24.529 

^  />'|2yoVdx/  j 


Ffyf 


J 


2 


can  be  regarded  as  the  nonlinear  stiffness  term  of  a  hard  spring  oscillator.  The 
forcing  term  for  the  first  mode  is: 


f(t)  = 


t)<fe(x)dx  . 


In  this  report,  we  shall  assume  f(t)  to  be  a  bandllmited  gaussiau  random  process 
of  white  noise  type  and  zero  average.  Of  greatest  interest  are  the  probability 
distributions  of  the  displacement  crests  and  extrema  of  a  hard  spring  oscillator 
under  gauss ian  forcing*,, 

1.  Displacement  crests 

The  average  number  of  positive-slope  crossings  of  the  level  y,  per  second,  is 
given  by  (ref.  5) 


pco 

n^(y)  =  /  uW(y,u)  du 

*'0 


(II-5) 


Writing  the  equations  of  motion  in  the  following  ^^system  form",  for  a  general 

g(y)  =  -  gC-y) : 


=  —  2au  —  g(y)  +  f 


(11-6) 


we  find  the  following  Fokker-Planck  (F-P)  equation  (ref.  6)  for  the  joint  distri¬ 
bution  of  y  and  u:  i.e.,  for  W(y,u,t) 


dw  ^  dw  ^  ,  I  n 


(II-?) 


assuming  white  noise  excitation  of  spectral  density  D. 

The  asymptotic  (i.e.,  stationary)  solution  of  this  F-P  equation  can  be  found 
by  analogy  with  a  quantum  mechanics  problem  (ref.  7)  to  be: 


W(y,u)  =Aexp[-  ^  (u2-t  gtOdC )]  = 
=  ^exp  [-  •j^  (u2+2/g(£)de)j 


(II-8) 


3 


where 


(II-9) 


i 

i 

I 

I 

i 

1 

i 


I 


D 

4a 


which  can  be  regarded  as  the  actual  mean  square  if  the  bandwidth  of  the  f  process 
is  sufficiently  large. 

The  expression  for  n*^(y)  now  becomes,  in  the  hard  spring  case: 


n'^(y)  =  A<r^e%p 


1/2  4  \ 

T*  > 


(II-IO) 


where  represents  the  average  number  of  positive-slope  zero  crossings, 


'y*  ,  and  can  be  expressed  as: 


V  *  -- 


[3.6256  -2>/2  l.2254,F,(f  .f  ;  z  )] 


(il-li) 


where 

7  = 

4/3  <r^  ~  ^  D 

after  making  a  change  of  variable  and  using  a  published  definite  integral  (ref.  8). 
This  integral  also  gives  alternate  expressions  for  the  moments  of  y,  which  are 
practical  for  larger  values  of  D,  by  interpolating  between  tabulated  values  for  the 
confluent  hyper  geometric  functions  in  Jahnke  and  Emde  (ref.  9). 

For  computing  the  relative  crest  distribution  we  need  not  know  ,  however. 

There  we  need:  ® 


c(y)  5 


n(v) 

n(0) 


n"(y) 

n"(0) 


(11-12) 


Thus,  if  displacements  exhibit  hard  spring  behavior  the  deviation  from  the  gaussian 
crest  distribution  (p  =  0)  should  be  discernible.  Experimental  displacement  crest 
measurements  do  not  indicate  appreciable  deviation  from  the  gaussian  crests  so  that, 
for  strain  calculations  at  least  (see  §IIB),  we  may  consider  the  y  process  as  ap- 
proximately  gaussian.  At  worst,  we  may  be  said  to  have  ’’equivalent- linearized”  the 
y  process.  More  incentive  to  use  such  an  approximation  scheme  will  be  provided  by 
the  fact  that  the  distribution  of  the  maxima  for  this  process  cannot  be  normalized 
in  terras  of  known  functions,  as  outlined  in  the  next  section. 

2.  Displacement  extrema 

The  probability  distribution  for  the  maxima  of  the  stationary  y  process  is 
given  by  (ref,  5): 


Wmaxly)  "  “  W(y,o,v)vdv 


(11-13) 
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where  W(y,u,v)  is  the  simultaneous  joint  probability  density  of  y,  u  =  ^  and 

V  =  ^  ;  and  A  .  is  the  normalization  constant  making W  (y)  dy  =  I, 

—  dt^  *  -  00 


If  f(t)  has  sufficiently  wide  bandwidth,  B,  we  may  treat  y  and  f,  and  u  and  f, 
as  approximately  statistically  independent,  since  y  and  u  will  fluctuate  much  more 
slowly  than  f  (here  we  have  tacitly  assumed  ergodicity  (ref.  10)  ).  Thus,  to  find 
V'maxCy)  we  perform  the  following  equivalent  integration  by  using  the  following 
result  obtained  from  Eq.  (II-6)  : 


W(y)  =-A,/w(y,0)W(f  )df 

max  V 

-oa 


=  v(f ) 


u«0 


u>0 


=  -“fy  f 


to  determine  the  limits  of  f.  We  use  the  W(y,u)  cited  earlier,  and  our  hypothesis 
implies  that : 


W(f  )  = 


I 


exp  - 


(ii-i^) 


We  let  a  £  =  DB  and  hence  assume  bandlimited  white  noise  for  a  sensible  calculation. 
Thus,  for  the  maxima: 


W  (y)=  A, 

max 


«*^|y +^y^  (i+erf  (  ‘vlyt.P-ii. 

2  \  \  y? 


20-2 
f 


(11-15) 


u 


The  above  formulation  can  be  extend'd  to  cover  the  arbitrary  odd-function  re¬ 
storing  force  of  formula  (11-6)  which  predicts  that 


W(y)  = 

max 


(y ) 

2  .2 


W(y,  0). 


(11-16) 


The  first  term  on  the  right  hand  side,,  i.e.  ,  2  S(y)  W(y,o)  can  be  interpreted 
as  the  distribution  of  the  excess  of  maxima  over  minima,  and  can  be  obtained  by  dif¬ 
ferentiating  the  crest  distribution .  Here  the  minima  are  treated  like  negative  maxima 
and  subtract  out  in  the  counting  process;  under  the  assumption  of  no  negative  minima, 
as  made  by  Lyon  (ref.  15),  this  can  be  interpreted  as  the  distribution  of  positive 
maxima.  The  other  two  terms  on  the  right  hand  side  in  Eq.  (11-16)  can  be  considered 
as  ’’correction”  terms  due  to  the  existence  of  positive  minima  and  also  negative 
maxima  of  displacement. 


The  y  process  is  statistically  symmetric  about  the  origin,  as  can  be  verified 
by  substituting  -y  for  y  in  the  equation  of  motion  for  y,  and  also  by  noting  that  the 
f  process  possesses  this  s^mimetry.  Thus  the  probability  density  of  the  minima  is 


5 


given  by: 


Wmin(y)-  Wmax('“y) 


Unfortunately,  difficulties  arise  when  it  is  attempted  to  find  the  normalization 
constant  Aj  for  the  maximum  distribution;  numerical  integration  appears  to  be 
necessary.  Apparently  only  the  linear  oscillator  case,  p  0,  yields  A(  in  closed 
form.  This,  plus  the  experimental  displacement  crest  measurements  already  cited, 
forces  us  to  approximate  the  y  process  by  a  gaussian.  Thus,  we  proceed  as  if  p=0, 
in  the  sequel. 

We  now  assume  that  the  joint  probability  density  of  y,  u,  and  v  is  of  the  form 
(ref  s  .  5,11)  : 


W(y,u.v)  =  -^r-==== 

(arr)  vf'(mofn4-  mlJnrig 


exp 


-  — f-yi  + 

rrvjV^+Zmgyv  + 

2  \  m2 

(11-17) 


which  is  the  gaussian  distribution  specialized  to  the  case  where  the  correlation 
matrix  is : 


<yy>  <yu>  <yv> 

1 - 

CVJ 

E 

1 

o 

o 

E 

L _ 

<uy>  <uu>  <uy> 

= 

0  rrig  0 

<vy>  <5/u>  <^vy> 

! 

1 

3 

!V 

O 

3 

1 _ 

as  demanded  by  the  stationarity  assumption.  Here 


For  this  case  the  normalized  variable  = 
bution  of  the  maxima  (ref.  11): 


possesses  the  following  distri- 


Wmax[77]  =  ^  «  exp  ^-  exp  ("V)  +  \') 

X 

where  the  parameter  t  can  be  found,  variously,  as: 


=  1  - 


rhe 


- 

V  *'mox  / 


=  1  -  (I  -2n 


(11-20) 


where  r  =  avg.no.  of  negative  maxima  per  second/y 


maj^ 


6 


Let  us  now  compare  the  above  result  for  a  general  gauss ian  process  with  the 

expression  for  _(y) , evaluated  at  6  =  0, 

^  max  ’ 

In  the  linear  oscillator  case: 


Stationarity  demands  that 
relations”  must  hold: 


-  2au  —  a|2y 


+  f 


(11-21) 


Thus,  the  following  "consistency 


<VU>  =  0 


-  2a 


<Ju> 


i  •  e  • , 


<fu>  =  2a  <ub> 

However,  from  formula  (II -9)  and  from  <r^  =  DB 


<tH> 


2  a 

\/4<<B 


(11-22) 


and  for  a  sufficiently  large  bandwidth  B,  f  and  u  may  be  regarded  as  uncorrelated; 
for  gaussian  statistics  this  means  statistical  independence.  For  p  =  0,  of  course 
(ref.  6), 


m. 


=  <i^=  0.2 


nri/ 


(11-23) 


Also 


^4  =<v^>=  4a2<:^u^>  +  4a  <uf>  + 


Since  u  and  f  may  be  regarded  as  practically  uncorrelated,  we  expect  that  y,  an  even 
more  regular  function  of  time  than  u,  is  uncorrelated  with  f .  Using  =  ci 

we  find  that; 

1714  =  DB  +  “'a^o  “  4a2m2 

which  may  be  rewritten  as. 


m4  =  ^4aa»|B  +  o.^  -  4a^wo^mo  .  (II-24) 
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T  hus , 


u> 


fm? 


+  4auiJ(B 


+40  ( B  -  a) 


U) 


2  -+  4a0 
0 


(11-25) 


for  sufficiently  large  B,  This  formula  for  e,  when  substituted  into  the  "exact** 
expression  for  W  ^(y)  evaluated  at  p  =  0  provides  the  proper  gauss ian  limit  for 
Wjnax^y^f  as  expc^cd. 

The  y  process,  then^  will  be  considered  to  be  gaussian  in  discussions  of  strain 
statistics  in  the  following  section* 


B*  Strain  equations 


The  bending  stress  at  the  bar  surface  is  given  by 


T,  = 


(11-26) 


where  h  is  the  half -thickness  of  the  bar*  For  the  fundamental  mode  of  a  clamped- 
clamped  bar: 


-Y„h  X,A 


coshX,x  +  cosx,x  —  0.9825  (sinhx,x  +  sinX,x)  y|(t).  (ii-27) 


Evaluated  at  the  center  of  the  bar  this  becomes 


while  at  3/8  X. 


27.197 


YohA 


^1 


the  experimental  location  of  the  strain  gages  : 


.777  T., 


(4) 


(11-28) 


(11-29) 


The  corresponding  ratio  of  the  corresponding  displacements,  for  the  fundamental  mode, 
is : 


y(i^-t)  =  .872y(-|-.t)  .  (11-30) 

The  tensile  stress  term,  constant  along  the  bar,  was  given  earlier  in  Eq.  (II-2). 

If  linear  stress -strain  relations  are  assumed,  then  the  total  strain  at  a  given 
point  on  the^urfacc  of  the  bar  is  proportional  to  the  sum  of  the  tensile  stress  and 
the  bending  stress  at  that  point,  and  can  be  written  in  the  following  manner: 

s  =  ay  +  by^  (ii-3i) 

where  y  is  a  convenient  displacement  variable. 
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Various  choices  of  displacement  variable  y  may  be  made;  y  may  represent  the 
displacement  at  x  =  V2  or  at  3/8j{  ,  with  the  necessary  changes  in  the  constants  a 
and  b.  We  may  then  write 


s{i>t)=  a,y(-|-)  + 


(11-32) 


with 


Q,  =  17.125-^ 


b,  =  2.446 


(11-33) 


A  theoretically  convenient  normalization  scheme  for  describing  the  statistical 
behavior  of  the  strain  involves  normalization  with  respect  to  the  rms  bend ing 
strain,  i.e,  ,  we  take  _  • 


Thus 


s(-^) 

a,<ry(f) 


^  +  y(4)’? 


1^2 


(11-34) 


where 


or 


>'li>  -=  =  .143=5^' 


(11-35) 


where 


y{¥)  -  - 

b, 


yi 

It)  1 

f  -X-  i 

V  2  / 

y\T^ 

^  8 

t777){.872)a, 


.7770, 

_  p  1 1  <^y('8 

h  ■  h 


(11-36) 


also. 


Al*)  - 


r(IL 

,677 


Experimentally,  of  course,  it  is  more  convenient  to  measure  the  rms  of  the  total 
strain,  or  rather  this  rms  value  without  inclusion  of  the  dc  component.  This  means 
that  a  strain  renormalization  process  is  necessary  before  the  theory  can  be  compared 
with  experiment;  a  discussion  of  this  appears  in  Appendix  I, 
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1.  strain  crest  statistics 


To  find  the  average  number  of  crossings  per  second  of  the  level  U  with  positive 
slope,  we  have  to  evaluate  the  analog  of  Eq.  (II- 5) 


where  £  =  However,  since  we  know  the  relation  between  o)  and  f)  ,  we  can  perform 

an  equivalent  ^  integration  instead* 


Since, 


(11-37) 


a  positive  slope  crossing  of  the  level  0  must  be  produced  by  either  of  two  possi¬ 
bilities  ,  either  i)  by **7  going  through  the  appropriate  value  above  -  ^  (to  be 
called  ‘^-1.(713)  later)  with  '7  >  0  or  ii)  by  'JJ  going  through  the  proper^Y  value  below 
-  i  (to  be  called  **7  (0))  and  7  <  0*  Thus: 

2y 

/OO  Q 

w('»fl[«’),u)udu  —  y^W('7lor),u)u  du  =  v^c(w)  .  (11-38) 


For  a  general  odd-function  restoring  force: 


n'^M  =  "o  |exp|j  +  exp  ~  j 


(11-39) 


For  a  hard-spring  oscillator  excited  by  white  noise: 


n‘*^(Tir)  =  *'o'sexp 


+  exp 


2‘^u^ 


2  cr 


v:hile  for  a  linear  oscillator  (p  =  O)  this  becomes: 


(11-40) 


=  "a  [- 


exp 


[- 


7^(®) 


(Il-i^l) 
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Thus,  the  avera^^.e  number  of  positive-slope  zero  strain  crossings  per  second, 
in  the  general  case:  becomes: 


n-^(o) 


exp 


2k 


g(e)dt 


(II-4?) 


where  ®’y/Y  is  a  constant  (i.e.,  is  independent  of  °’y)  we  recall.  For  a  linear 
oscillator : 


n'^(o)  =  i/t 


(11-43) 


2.  Ebctrema  of  strain 

At  the  extrema  of  strain,  from  Eq.  (11-37) 


0  =  •*■  =  (1  +  Zyr])'n 

Thus,  either  '>7=0  (displacement  has  an  extremum)  or  V  =  ~ latter  produce 
the  absolute  minimum  value  of  w  ,  =  -  1/4y,  To  investigaTr'e  the  specific  nature 

of  these  extrema  we  need  the  second  derivative; 


(  I  +  2YV) 


d^V 

dt2 


(Il-44> 


Thus,  extrema  of  V  occurring  for  t7  above  -  l/^Y  (i.e.,  extrema  of  y  above 
the  appropriate  constant  multiple,  of  h)  produce  extrema  of  having  the  same  sign 

of  '9  ;i.e.,  maxima  of  tj  go  over  into  maxima  of  w  ,  and  minima  of  17  go  over  into 

minima  of  w  .  If  y  is  small  enough  this  is  all  we  need  consider;  this  is  the 
approach  of  ref.  4,  18. 

However,  the  sign  of  the  second  derivative  of  tx>  is  opposite  that  of  V  for 
extrema  of  ro  occurring  for  ^7  below  -  l/^y;  maxima  of  go  into  minima  of  txr  ,  and 
minima  of  17  go  into  maxima  of  iir  in  this  range. 

Extrema  of  t]  occurring  precisely  at  '>7  =  -  1/2y  produce  an  inflection  point 
of  Tzr  since  bothiaT  =  *ar  =  0,  while  '*7  crossing  of  the  level  V  =  -  1/2y  (making 
ar  =  0)  with  “i?  ^  0  produce  minima  of  vr  ,  since -or  =  2y(  V  ^  0. 

Thus,  we  are  able  to  designate  extrema  of -or  as  appropriate  extrema  of  17  or  ^ 
as  crossings  of  the  level  V  =  -  1/2y.  We  may  then  find  the  statistical  properties 
of  the  extrema  of  nr  from  knowledge  of  the  properties  of  the  17  process. 

a)  The  maxima  of  Tir  are  thus  caused  either  by  maxima  of  V  above  -  l/2Yf  or 
by  minima  of  V  occurring  below  -  1/2y.  In  particular  the  number  of  maxima  of  ■®‘ 
per  second  is  given  by  the  sum  of  the  average  number  of  maxima  of  V  above  -  1/2y 
and  the  average  number  of  minima  of  ^  below  -  l/2x*  ^ov  gaussian  V  this  is: 


V 


*'mox 

V 


max 

■or 


(11-45) 


i 


The  equation 


r  V  ^  +  V 


=  0  admits  to  two  roots  : 


I 

( 

1 

i 

I 

i 

{ 

I 

i 


17--^ 

-Jr  - 


(11-46) 


Evidently,  17  +  is  in  the  range  of  r)  values  above  -  1/2y,  while  71  is  always  below 
-  1/2y,  as  anticipated  earlier. 

The  probability  distribution  for  the  maxima  of  ur  can  be  obtained  by  performing 
the  appropriate  *>7  integrations  once  77  )  and  Wj^in(*'9)  known.  From: 


) 


(11-47) 


we  deduce,  after  changing  the  integration  variable  to  w  ,  that: 


=  a' 


w, 


minb'H  + 


v/TTTyi 


(11-48) 


where  hf  is  the  normalization  constant.  Here  ^  ^  “  ^rain^"^^  that 


/l  +  4XW- 


(11-49) 


The  second  terra  on  the  right  hand  side  of  Eq,  (11-49)  is  the  only  term  of  importance 
for  small  y>  ^nd  hence  is  the  only  one  appearing  in  Lyon's  analysis  (ref.  4),  His 
analysis  appears  formally  different  because  he  normalizes  the  y  process  with  respect 
to  its  ’’linear"  rra.s  value.  The  analysis  of  P,  W.  Smith,  Jr.  (ref.  17)  also  assumes 
small  Y  but  uses,  in  effect,  the  previously  discussed  first  term  of  Eq,  (11-15)  as 
W 

max 

b)  Similarly,  the  minima  of  tir  can  be  attributed  to  either  minima  of  17 
occurring  for  77  >  -  1/2y»  or  to  ri\axima  of  77  ,  occurring  for  V  <  ~  i/^Yt 
crossings  of  the  level  -  1/2y  by  the  '7  process. 

The  contribution  of  the  extrema  of  ^*7  to  the  distribution  of  the  minima  of  txt  can 
be  found  as 


^3) 


=  A"| 


ilTlnr) 

W, 


T7-(nr) 


max 


4- 


Kin 


(11-50) 
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which  leads  to: 


=  A" 


I  4-  4yw 


(11-51) 


where  A"  is  the  normalization  constant  making 


^minWd®  =  I  -P(m^,=  -  '  ) 


In  this  case  ''’njijjC*?  )  =  '‘^tnax^  ^ 


‘\ACi„(a7)  =  A" 


«/  W„ 


v/1  +  4/  W 


(11-52) 


Note  again  that  only  the  second  term  of  the  right  hand  side  of  Eq.  (11-52)  is 
significant  when  the  nonlinearity  parameter  y  is  small  (ref,  4) • 

The  probability  that  a  minimum  of -or  is  caused  by  a  crossing  of  the  level 
-  1/2y  by  17  can  be  found  simply  in  terms  of  the  average  number  of  minima  of  or  per 
second  caused  by  extrema  of,  i.e.  the  number  of  crossings  of  the  level 

1  -  •  l/2r  per  second,  i,c.,  n(-  1/2y)  =  2n+(-  1/2y) 


in  27  ^ 


(11-53) 


For  the  gaussian-derived  process; 


=  ''mox 


(11-54) 


with  I/- 


Also : 


<-w)-  2*';exp(--^ 


(11-55) 


P(®^min  =  “  )  = 


'J  \  ~ 


(11-56) 
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Notei  «8  Y  0  lain)  ^  Of  while  for  y  a>  it  becomes 


Using  this  result,  we  find  that 


+  I 


^  I 


I  +y,-  «2  'exp(- 

The  distribution  of  the  minima  obviously  has  a  jump  of  magnitude  P(‘®’njin  ”  *  i/^Y) 
at  =  -  1/4y;  and,  of  course,  there  are  no  minima  of  ur  below  -  1/4y* 

For  experimental  verification,  of  particular  interest  are  the  distributions 
for  negative  minima  and  for  positive  maxima;  i.e*,  we  want: 


and 


P(«g«nox;»  ^  ) 
P  (wmox^  O  ) 


(11-57) 


P_(wmln  ^ 


P  (”min  ^  ^ 

P  ( <  0  ) 


(11-58) 


Fpr  the  gaussian-derived  square  law  process  (ref.  12): 


P  (®mox  ^ 


2  -h 


e*p(-  “iV^) 
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Of  course.  >  -  IAy)  =  1.  From  this  ,  we  get  that : 

inax  ^ 

P|.(»«0K  ?  w)  -  +  /)  -  <2  - [*  ■*■  erf^ 

-  ‘'ir. 

*  ATT^ 


(11-60) 


Simplifications  can  be  performed  on  these  expressions  for  certain  cases:  If 
Y  is  small  (i.e. ,  y  «  1) , 

P+{'^mo*^  ®)  ■  ~  J{  -  exp^-  ^  [l  + 

(11-61) 


while  for  e  -  0, 

P+(«rtnox  >  O) 


1  - 


(11-62) 


Nov;  for  small  e  (i.e.,  e  «  1)  we  have 


P4-.,.sw)  •  ^  A^  eKp(-  [l  +  erf(:^' 


2A^e.p(-^)- 


(i  /rrT?  2xp(--4.)  X-  A^T^) 
V  2/^  '  / 


(11-63) 


Also,  after  straightforward  but  tedious  integrations  we  find  that  (ref.  12)  ; 

P(«^min  w) 


2v/i  -  'exp(--^ 


^^«p( -%>-)[,  -  e„(^  ^)j 


(11-64) 


’?_(»)  yi-  «2 


2  -*-  2  A-  «  ^  exp - '—p 
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As  expected,  PC'®’ ^  -  IAy)  =  0  and  ^  +  co)  =  1.  Using  this,  we  find 

that : 

R(«ini»<">-)  =  I  a/P?"  eKp(--j^)  <■ 

+  exp(--^)[l  -  )]-.  /ir^ exp(- %-' )  X 

X  [|  +  |2yPT^exp(-^-^)^  +  ert  (^^) 

+  '/^«p(-A)['  - 

c<  1,  this  simplifies  to: 

>  =  '  -  ^,xp(-^fi)[l 

(^’^)l)/ 


For  small  Yt  T  i»  this  simplifies  to: 

P-(®min  ^ 


erf 


2'  €  yjj 

When  both  y  e  are  small  compared  with  unity, 


^  I 


P-{®min$^)  = 


But,  if  e  =  0  and  y  is  arbitrary;  then 


- ^  /  ^^(jBr)\ 

Z/A  exp^-  2  I 


exi 


-I-  I 


P(-^)  4  e.p(-#^) 


ex 


p(--^)+  ' 


(11-66) 


(IT-67) 


(11-68) 


This  result  differs  from  the  crest  distribution  for  large  y.  However,  for  larger  y 
this  expression  becomes  of  academic  interest  because  of  the  proportionately  few 
minima  occurring  for  negative  ur  •  The  P.  curves  depend  more  critically  on  y  than  do 
the  curves,  owing  to  the  presence  of  a  cutoff  for  ^  at  nr  =  -  l/4y. 

Limitations  in  the  theory 

The  range  of  validity  of  the  single-mode  model  employed  here  is  open  to  debate. 
It  is  valid  only  if  the  shape  of  the  eigen-function  persists  for  large  amplitudes^ 
and  if  contributions  to  the  motion  from  other  modes  can  be  disregarded.  In  addition, 
numerical  integration  is  necessary  to  make  use  of  the  more  exact  expression  for 
W^x^y)  rather  than  gauss ian  approximations. 
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Unfortunately,  it  is  known  from  experiment  (refs,  13,1^)  that  mode  coupLin;^ 
cannot  be  neglected  at  larger  amplitudes.  Also  the  experimental  values  of  the  non¬ 
linear  stiffness  differ  drastically  from  the  theoretical  values  (ref.  14),  Moreover, 
the  mode  shape  changes  at  larger  amplitudes  and  the  modal  damping  must  be  regarded  as 
nonlinear.  In  short,  this  analysis  breaks  down  for  large  y-  However,  it  is  expected 
that  analysis  here  developed  will  have  validity  for  larger  y  values  than  preceding 
ones  (refs.  4,17)  because  of  the  inclusion  of  strain  maxima  (or  minima)  that  are  not 
caused  by  displacement  maxima  (or  minima,  respectively). 
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Section  III.  Configuration  and  Instrumentation 


A.  Design  Considerations 


In  order  to  measure  effects  of  dynamical  nonlinearity  on  stress  statistics  it 
was  necessary  to  design  a  clamped- clamped  bar  and  driver  system  such  that  the  tensile 
stress  could  be  of  comparable  magnitude  to  the  bending  stress  at  the  surface  of  the 
bar.  The  bending  stress  at  the  surface  is  as  given  in  £q.  (11-26), 


-  Yoh 


where  h  is  the  half  thickness  of  the  bar. 

The  displacement  of  the  bar  as  a  function  of  x  is  of  the  form 


=  A(coshX„x  -  cosXnx)  -t  B(sinhX„x  -  sinX^x) 


(iii-i) 


where  A  and  B  are  constants,  n  is  the  mode  number,  and 


X^ 


For  the  fundamental  mode,  Eq.  (III-l)  becomes 

^,(x)  =  A  |coshX,x  -  cosX,x  —  0.9825  (sinhx,x  -  sinX,x)J 


X.  =  1.004 


iir 


~2T 

=  -f 


with  A  chosen  so  thot 


i.e.  A  =  .70639 


(III-2) 


By  substituting  Eq,  (III-2)  into  (II-2)  and  (11-26),  the  following  stress  relations 
are  obtained: 


6.l703YoA"  2,^, 

T,  =  - -  y,(t) 


(III-3) 


and 


Tg  =  -  X,Yj,hA  j^coshX,x  -»■  cosX,x  —  0.9825{sinhX,x  ■+  sinX,x)jy,(t)  -  (m- 


4) 
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The  bending  stress  evaluated  at  the  center  of  the  bar  becomes  Eq*  (11-28) 


Tg  =  27.1972  y(t) 


and  the  displacement  at  the  center  of  the  bar  is 


=  I.5882A 


Eliminating  A, Eqs.  (II I -3)  and  (III-4)  become 


3.8851  Y„  ..  ,  H 


A 


y  (  ^  .t) 


and 


l7J245X>h  n 

nZ  y V  2  » ’ ^ 


(III-5) 


(III-6) 


(III-7) 


Thus,  the  condition  for  these  stresses  to  be  of  comparable  magnitude  is  that  the 
amplitude  of  vibration  at  the  center  of  the  bar  be  of  the  order  of  the  ha  If -thickness 
of  th  e  bar . 

B.  Clamped-Clamped  Bar 

The  clamped- clamped  bar,  sho\^  in  Fig.  1,  was  machined  from  a  solid  piece  of 
2024-T4  aluminum  so  that  the  bar  and  its  boundaries  were  integral  one  with  the  other. 
This  configuration  reduces  clamping  losses  and  hopefully  effects  of  unwanted  non- 
linearities  at  the  boundaries.  The  bar  was  excited  by  driving  an  aluminum  voice 
coil,  mounted  directly  on  the  bar,  with  a  field  provided  by  a  dynamic  speaker  field 
coil.  The  force  thus  obtained  was  proportional  to  the  current  through  the  voice 
coil  over  the  current  range  used.  This  driver  unit  was  capable  of  exerting  a  force 
of  up  to  5  newtons  per  ampere.  The  stress  was  measured  with  strain  gages  in  the 
circuit  shown  in  Fig.  2.  This  circuit  makes  possible  the  separation  of  the  two 
types  of  stress  in  addition  to  determining  their  net  effect.  Such  procedures  arc 
accomplished  by  mounting  the  gages  on  opposite  surfaces  of  the  bar  and  then  adding 
or  subtracting  their  outputs  as  illustrated  in  Fig.  2.  The  strain  gages  were  placed 
at  a  distance  of  3JI/8  from  one  end.  The  rms  displacement  was  measured  with  a 
dis  placement  probe. 

The  response  of  the  bar  under  sinusoidal  excitation  is  shown  in  Fig.  3.  The 
bar  was  excited  at  its  fundamental  frequency,  135  cps.  Figure  3a  illustrates  the 
tensile  stress.  Fig.  3b  the  bending  stress,  and  Fig.  3c  and  Fig.  3d  are  the  net 
effects.  It  may  be  noted  that  there  is  a  rectification  of  minima,  hence  a  tendency 
to  reduce  the  amplitude  of  stress  minima  while  the  amplitude  of  stress  maxima  is 
augmented.  The  appendix  gives  the  modifications  introduced  by  the  removal  of  dc  strain. 

The  block  diagram  in  Fig.  4  illustrates  the  method  of  measuring  the  distribution 
of  crests.  The  random  forcing  was  obtained  from  the  noise  generated  in  a  photo¬ 
multiplier  tube  which  was  amplified  sufficiently  to  drive  a  power  amplifier  which  in 
turn  excited  the  electro-magnetic  transducer.  A  filter  in  the  amplifier  chain 
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shaped  the  excitation  spectrum  such  that  only  the  fundamental  mode  of  the  bar  was 
excited.  The  response  of  the  bar  was  detected  with  the  strain  gages.  The  output 
of  the  strain-gage  circuit  was  amplified  and  fed  into  the  amplitude  discriminator 
which  produced  a  large  pulse  when  a  preset  amplitude  was  exceeded  with  positive 
slope.  This  preset  level  was  set  with  a  test  oscillator,  of  known  output,  in  the 
circuit.  The  number  of  crossings  of  this  level  was  then  counted  on  an  electronic 
digital  counter.  This  procedure  was  repeated  for  enough  levels  to  determine  ex¬ 
perimentally  a  sample  frequency  distribution  of  the  positive  crests.  The  distri¬ 
bution  of  the  negative  crests  was  found  by  inverting  the  signal  into  the discrimin^or 
and  repeating  the  entire  process. 

C.  Instrumentation 

Many  of  the  measuring  instruments  used  in  these  experiments  were  standard 
commercial  units  as  shown  in  the  block  diagram  of  Fig.  4,  The  differentiator  and 
amplitude  discriminator  circuits  were  developed  from  data  furnished  hy  the 
G.  A,  Philbrick  Co.  and  are  shown  in  Fig.  5.  Square  wave  tests  of  the  differentiator 
revealed  good  performance  well  into  the  upper  audio  range  where  its  action  was  pur¬ 
posely  degraded  to  reduce  noise  problems.  Differentiated  noise  from  the  amplifiers 
masked  the  desired  random  response  signals  unless  some  bypassing  was  done.  The 
discriminator  circuit  is  derived  from  a  Philbrick  operational  amplifier  with  one  of 
the  differential  inputs  biased  by  a  three  volt  battery.  When  the  potential  of  the 
other  input  exceeds  three  volts,  the  discriminator  output  goes  to  70  volts.  The 
device  operates  at  frequencies  up  to  100  kc  and  its  output  is  recorded  by  the  digital 
counter.  Any  significant  error  introduced  by  the  discriminator  circuit  is  due  to  its 
50  mv  threshold  which  represents  a  small  correction  appreciable  only  at  low  preset 
levels . 

The  rms  level  of  the  strain  voltage  was  measured  by  using  a  capacitor  bank  of 
30,000  mfd  to  integrate  the  signal  feeding  the  Ballantine  320  Voltmeter  which 
operates  as  a  square- law  detector.  The  linearity  of  the  amplifying  chains  used  was 
excellent  for  the  signal  levels  used.  Noise  pickup  by  the  strain  gage  wiring  was 
reduced  by  electric  shielding. 
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Section  IV*  Results  and  Conclusions 


The  theoretical  distribution  of  positive  maxima  is  compared  to  the  theoretical 
distribution  of  positive  crests  in  Fig.  6  for  the  nonlinearity  parameter  0*0475. 
The  Rayleigh  distribution  which  corresponds  to  y  =  0  is  also  indicated.  For  e  =  0 
the  maxima  distribution  coincides  with  the  crest  distribution.  As  e  increases  the 
maxima  distribution  falls  below  the  crest  distribution.  Similarly,  in  Fig.  7,  the 
negative  minima  and  the  negative  crest  distributions  coincide  when  e  =  0  while  the 
minima  distribution  falls  belov;  the  negative  crest  distribution  as  e  increases. 

The  theoretical  and  experimental  distributions  of  positive  crests  for  two  values 
of  the  nonlinearity  parameter,  namely,  y  =  0  and  y  =  0.082  are  shown  in  Fig.  8,  and 
the  distributions  for  y  =  0,113  is  shown  in  Fig.  9.  The  value  y  =  0  corresponds  to 
linear  motion  for  which  the  crest  distribution  is  Rayleigh.  The  distributions  re¬ 
sulting  from  nonlinear  motion  clearly  deviate  from  the  Rayleigh  distribution.  As 
the  nonlinearity  parameter  y  increases  the  deviation  from  the  Rayleigh  distribution 
increases.  At  small  excitation  levels,  where  the  tensile  stress  is  very  much  less 
than  the  bending  stress,  the  experimental  points  follow  the  Rayleigh  curve  within 
the  limits  of  experimental  accuracy  while  at  higher  excitation  levels  the  experimental 
points  deviate  somewhat  from  the  theoretical  curv*es  particularly  at  small  values  of 
07  •  At  high  values  of  jjt  ,  however,  both  theory  and  experiment  indicate  that  the 
positive  crest  distribution  falls  above  that  expected  for  linear  motion.  The  number 
of  positive  crests  occurring  above  two  times  the  rms  stress  is  appreciably  above  the 
number  occurring  in  the  same  region  for  the  linear  case. 

The  theoretical  and  experimental  distributions  of  negative  crests  for  y  =  0 
and  y  =  0.082  are  shown  in  Fig.  10  while  the  distribution  for  y  =  0.113  is  shown  in 
Fig.  11.  The  Rayleigh  distribution  is  again  shown  for  comparison  purposes.  The 
distributions  for  nonlinear  motion  fall  below  that  for  linear  motion,  particularly 
for  stress  levels  above  the  rms  stress.  Hence,  there  are  fewer  negative  crests  of 
large  amplitude  in  the  nonlinear  case.  It  therefore  may  be  tentatively  concluded 
that  the  Rayleigh  distribution  predicts  too  few  positive  crests  and  too  many  negative 
crests  for  high  stress  levels  under  conditions  where  dynamical  nonlinearities  of 
this  kind  are  important. 
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Appendix  I.  Normalization  of  Strain 

The  expression  for  strain  used  in  this  report  was  Eq.  (11-31) 

s  *  ay  by^  . 

Thus,  the  dc  component  of  strain  is  given  by 


<s>  =  b<y2>  i  b<r/  .  (A-i) 

EJxperimentally,  the  dc  component  is  often  removed  from  the  strain  (as  was  done  here) 
and  the  rms  value  of  the  remaining  signal  is  measured.  In  other  words,  the  experi¬ 
mentally  measured  strain  is 


sis  —  -  ay  by^  —  b<ry^ 

and  the  experimental  rms  strain  value  is  : 


(A-?) 


— 


4  2  2 

In  particular,  for  a  gauss ian  y  process  ^yj>  =3  ^ 


(A.3) 


and 


Srmt  =  a<r„  /l  +  2y2 


(A-4) 


where,  as  before 


For  a  hard  spring  oscillator, 


however  (ref. 


17) 


<V^)  •  ^(<y%  -  <v‘>)  <*-” 


which  can  be  rewritten  as  (ref.  13) 

<y^  =  4z  (<^5>o  - 

=  4z  [l 


<y\  “■y  ) 


-  4/T 


1.2254, Fi  ( 


3.6256 


i.etze/r  ,F- (y;j^z)  [- 

2 .4508 /r  .  F?  V  z  )  I  / 


It  thus  appears  preferable,  when  applying  the  theory,  to  renormalize 
strain  data  to  be  consistent  with  the  theoretical  conditions,  rather 
more  ^’realistic*'  normalization  scheme  in  the  theory. 


the  experimental 
than  to  use  a 
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The  experimentally-measured  strain  has  the  minimum  value: 


^min 


4b 


(I  + 


-  :^-(i  t  4/^) 
4b 


(A.6) 


Thus  the  experimental  strain,  unlike  the  total  strain,  does  not  have  an  absolute 
minimum  value  at  -  a^/^t);  rather  its  minimum  value  increases  negatively  as  y  increases. 


Using  the  normalization  scheme  of  the  theory  in  this  report,  the  normalized 
experimental  strain  becomes : 

^  YV^  -  Y 

where,  as  before,  7)  =  y/cry.  Thus  w  is  shifted  down  by  the  amount  y»  and  the 
minimum  value  is 


"min 


I 

4r 


(  I  +  4/^) 


(A-8) 


Hov;ever,  using  the  gaussian-type  normalization  scheme  of  Eq,  (A-4)  for  the 
experiment a I  strain 


+  zr 


and 


/ 


y  I  +  ZY^'  J I  +  ZY^' 


- 


mm 


I  ( I  +  4y  ) 


4y 


J\  +  ZY^ 


(A.9) 


(A-IO) 


Y  — ^  00,  -4  min  =  -  .707;  thus  for  large  displacements  ^  ^ixi  experi¬ 
mentally  will  occur  at  approximately  -  .7S  ,  and  x^ill  become  more  negative  with, 

increasing  rms  strain  unlike  the  "total"  strain  which  has  an  absolute  minimum  at 
S  =  -  av4b  always,  provided  that  our  single-mode  model  is  valid. 


This  shift  in  the  normalized  strain  to  with  dc  removed,  as  contrasted  v/ith  or  , 
should  be  accounted  for  when  using  the  formulas  of  this  report  at  larger  y  values. 
This  can  be  accomplished  easily  since 


-or  =  w  -t-  X 


(A- 11) 


obviously.  This  enables  us  to  change  or  reinterpret  all  the  strain  formulas  in  this 
report  straightforwardly.  For  example,  the  previously  designat ed  zero  crossing 
forraulas  (for  m  )  nov?  give  the  crossings  of  the  level  3j  =  -  y,  ^ 

becomes  ^  -  y)  ,  etc. 

lUoX 
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For  the  Y  range  of  expcrinental  interest  in  this  report  such  corrections  arc- 
insignif  icant  and  we  can  use  ibt  and  #  interchanr^eably ,  Si.i:nilarl37’  here  we  can  inter¬ 
change  the  theoretical  and  the  experimental  rms  values  of  strain  in  the  norrnali- 
zation  process  for  strain,  i.e#,  we  can  interchange  &  and  ^  .  As  long  as  we  deal 
with  ^aussian  y  processes  and  os  long  as  our  mathematical  model  is  valid  the 
modifications  entering  when  y  is  larger  are  straightf orward  (ref.  12). 
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Coefficient  giving  bending  strain  as  ay 
Constant  =  0,7065 

Coefficient  giving  membrane  strain  as 
Bandxvidth 
Spectral  density 
f (t)  Normalized  effective  force 
h  Half  thickness  of  bar 
Length  of  bar 
Mass  of  bar  =  pSL 

Rate  of  exceedance  of  positive  maxima 

Rate  of  exceedance  of  negative  minima 
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X 

M 


Strain  =  ay  +  by^ 

Cross  sectional  area  of  bar 
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d\ 
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d^y  _  ^ 

dt  dt 


(y)  Probability  density  of  the  maxima  of  y 


W 

max 

W^in(y)  Probability  density  of  the  minima  of  y 
W(y,u)  Joint  simultaneous  prob,  density  of  u  and  u 

W(y,u,v)  Joint  simultaneous  prob,  density  of  y,  u,  and  v 

Probability  density  of  maxima  of  normalized  strain 
^^min^^^  Probability  density  of  minima  of  normalized  strain 
y  Displacement  of  bar 

Nonlinearity  coefficient 
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Y 

V 

e 


Nonlinearity  parameter  =  — 

.  .  y 

Normalized  displacement  variable  = 

Normalized  rms  width  of  power  spectrum 
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Pe 

<r 
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1 

'/l  +  4y^ 
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1 

\ 

«/l  +  4yC5 

K 

Radius  of 

gyration 

Wave,  number  of  fundamental  mode 
Average  number  of  maxima  of  13  per  second 

Average  number  of  positive- slope  zero  crossings  of  y(or  ^  )  per  seconds 

Mass  per  unit  length 

RMS  displacement  of  nonlinear  bar 

Normalized  spatial  displacement  eigenfunction  of  first  mode 
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Normalized  strain  variable  = 
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Fig.  1  Diagram  of  Clamped-Clamped  Bar 


Fig.  2  Strain  Gage  Circuit  for  Response  Measurements 


29 


v/\A^\/\/\ 

(a) 


(C) 


Fig.  3  Sinusoidal  Stress  Response  of  Bar 


When  ein>3v,  6Qy|  =  70v 

K2~W  is  a  Philbrick  differential  amplifier 

Discriminator 


0.001 


Differentiator 


Fig.  5  Amplitude  Discriminator  and  Differentiator 

Circuits  Employed  in  Probability  Distribution 
Measurements 
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Positive  Maxima  and  Positive  Crest  Distributions  for  y  =  0,0475 


Crests  and  Minima 


Negative  Minima  and  Negative  Crest  Distributions  for  y  =  0.0475 


Fig.  10  Negative  Crest  Distribution  for  Y  =  Y  =  0.082 
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However,  when  the  excitation  level  of  the 
I  clamped -clamped  continuous  bar  is  sufficiently 
high  so  that  the  tensile  strain  becomes  com- 
I  parable  with  the  bending  strain,  then  the 
I  strain  crest  distribution  no  longer  follows 
the  Rayleigh  prediction*  At  high  strain 
I  levels  the  distribution  of  positive  crests  as 
I  V7ell  as  maxima  is  greater  than  the  Rayleigh 
I  prediction  and  the  distribution  of  negative 
crests  as  well  as  minima  is  less.  The  dis¬ 
tribution  of  positive  maxima  falls  below  the 
I  positive  crest  distribution  as  the  Q  of  the 
I  system  decreases.  Similarly  the  distribution 
of  negative  minima  falls  below  the  negative 
I  crest  distribution  as  the  Q  decreases. 
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